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1 Introduction 

The discrete Fourier transform (DFT) and algorithms for its efficient compu- 
tation (FFT) enjoy an enormous range of applications. One might roughly 
divide such applications into the analysis of data which is sampled in time 
or in space. Applications involving spatially sampled data are basic in nu- 
merical analysis, since constant coefficient partial differential equations and 
some of their discrete analogs, including the heat, wave, Schrodinger, and 
beam equations in one space dimension, have exact solutions in Fourier se- 
ries. Other applications of the DFT involving the analysis of spatially sam- 
pled data include noise removal, data compression, data interpolation, and 
approximation of functions. These applications often take advantage of the 
tight and explicit linkage binding the harmonic analysis of uniformly sampled 
data with the harmonic analysis of periodic functions defined on M. 

Many aspects of Fourier analysis can be developed through the spec- 
tral theory of the second derivative operator —D 2 = —d 2 /dx 2 acting on the 
Hilbert space of square integrable functions with period 1. An orthogonal 
basis of eigenfunctions is given by exp(27rikx) . If N uniformly spaced sam- 
ples x n = n/N for n = 0, . . . , N — 1 are considered, the restriction of the 
exponential functions to the sample points 



again provides an orthogonal basis for functions on the set {x n }. These 
sampled exponential functions are also eigenfunctions for 'local' operators 
such as A^r, which acts by 



index arithmetic being carried out modulo N. The fact that efficient FFT 
algorithms are available for sampled data is not a consequence of abstract 
spectral theory, but relies on the exponential form of the eigenfunctions and 
the arithmetic progression of the frequencies k. 

Our discussion now shifts to graphs and their refinements, with the circle 
and its uniform sampling serving as an example. Discrete graphs with Ny < 
oo vertices and Ns edges have a well developed spectral theory [7] based on 
the adjacency or Laplace operators acting on functions defined on the vertex 
set V. While these ideas play a role here, the lead actor is a continuous graph, 



exp(2mkn/N), k, n = 0, . . . , N - 1, 
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also known as a metric graph or, when a differential operator is emphasized, 
a quantum graph. Here the edges of a graph are identified with intervals 
[a n , b n ] C R, functions are elements of the Hilbert space Q) n L 2 [a n , b n ], and 
differential operators such as — D 2 provide a basis for both physical modeling 
and harmonic analysis. (The terminology for such 'continuous graphs' is not 
settled, with the terms topological graphs or networks also used.) 

The natural and technological worlds offer numerous opportunities for 
graph modeling, where applications such as those mentioned above can be 
considered with data sampled from continuous graphs. These include the 
sampling of populations along river systems, and the description of traf- 
fic density on road networks. Biological systems transport nutrients, waste, 
heat, and pressure waves through vascular networks, and electrochemical sig- 
nals through natural neural networks. Elaborate networks are manufactured 
in microelectronics, and for microfluidic laboratories on chips. 

Despite a rapidly growing literature on 'quantum graphs' [U [H [9], and 
some work [2J generalizing the classical theory of Fourier series to functions 
defined on a continuous graph, there is a very limited literature considering 
the harmonic analysis of sampled data for these geometric objects. The 
papers [T3l [H] consider sampling on continuous graphs with virtually no 
length restrictions on the edges. In contrast, this work exploits the features 
appearing when graphs have equal length edges, and obtains the following 
conclusion: every finite continuous graph with edges of equal length admits a 
family of DFTs closely analogous to that of the circle, and an FFT algorithm 
for their efficient computation. 

To unify the presentation a bit, we focus on simple discrete graphs, with- 
out loops or multiple edges between vertices, whose vertices have degree at 
least two. Given a continuous version of a nonsimple graph without bound- 
ary vertices, one can insert additional vertices of degree 2 to reduce to the 
simple case. This insertion of 'invisible' vertices can also be used to reduce 
graphs whose edge lengths are integer multiples of a common value to the 
equal length case. Algorithmically, this reduction increases the size of the 
discrete graph spectral problem whose solution is an important component 
of the DFT. 

There are four subsequent sections in the paper. The second section starts 
with a review of differential operators on continuous graphs. This leads to 
a more detailed discussion of the spectral theory of the standard Laplace 
differential operator on continuous graphs with equal length edges. Most of 
the material in this section was previously known [SJ [TT] , although Theorem 
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2.9 appears to be new. The third section explores the linked spectral theory 
of Laplacians for continuous graphs and their uniformly sampled subgraphs. 
The fourth section shows that efficient algorithms for Fourier analysis are 
available. The final section presents a simple example where many of the 
computations can be done 'by hand'. 



2 Spectra of continuous graph Laplacians 

2.1 Continuous graphs and the standard Laplacian 

In this work a discrete graph Q will be finite and simple, with a vertexset V 
having N v points, and a set S of N £ edges. Each edge e has a positive length 
l e . It is convenient to number the edges. For n = 1, . . . , Ng the edge e n is 
identified with a real interval [a n , b n ] of length /„. The resulting topological 
graph is also denoted Q. In an obvious fashion one may extend the standard 
metric and Lebesgue measure from the edges to Q. 

The identification of graph edges and intervals allows us to define the 
Hilbert space 

with the inner product 

(f,g) = / f~9 = ^2 / fni X )9n{x) dx, f n \ [d n ,b n ] -> C. 

We will employ the standard vertex conditions, requiring continuity at the 
vertices, 

lim f(x)= lim f(x), e(i),e(j)~v, (2.1) 
and the derivative condition 

d "fn(v) = 0. (2.2) 

Here the derivative d u f n (v) is f' n {a n ), respectively —f' n {b n ) if a n , respectively 
b n , is identified with v. 

The standard vertex conditions are used to define the Laplacian Cf = 
—f" of the continuous graph Q. Let V max denote the set of functions / G 
L 2 (Q) with /' absolutely continuous on each e n , and /" G L 2 (Q). The domain 
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V of C is then the set of / G V max satisfying the standard conditions (12.11) 
and (12.21) . From the classical theory of ordinary differential operators [T2"l p. 
S123] one knows that L is self adjoint with compact resolvent. This operator 
is nonnegative, with being a simple eigenvalue if Q is connected. Writing 
eigenvalues with multiplicity, the spectrum is thus a sequence = A < Ai < 
A2 < . . • , and there is a orthonormal basis of eigenfunctions. 

We now impose the additional requirement that all edges of Q have equal 
length. In this situation the operator £ will be denoted Aoo. Initially the 
edge lengths are taken to be 1, but the rescaling £ = x/L converts the system 
D 2 X Y = XY to 

D\Y = fiY, < £ < 1, n = L 2 \. (2.3) 

Values of eigenfunctions at graph vertices are not effected, and derivatives 
are scaled by L. Thus eigenfunctions with eigenvalue A for the standard 
Laplacian on the graph with edge lengths L are taken to eigenfunctions with 
eigenvalue L 2 X for the standard Laplacian on the graph with edge lengths 1. 

The spectrum of A^ is tightly linked to the discrete graph Q and the 
spectrum of the discrete Laplacian of combinatorial graph theory. Let E(X) 
denote the eigenspace of A^ for the eigenvalue A. The connections between 
E(X) and the discrete graph Q arise in two ways, depending on whether or 
not A G {n 2 7r 2 I n = 1, 2, 3, . . . }. Before exploring this dichotomy, we note 
the following result. 

Proposition 2.1. Ifu>0, 

dim E(lu 2 ) = dim E([u + 2wr] 2 ), n = 1,2,3, ... . 

Proof. Suppose n > 1 and Y is an eigenfunction with eigenvalue uj 2 . On each 
edge 

y(x) = A cos(ux) + B sin(u;x). 

The function 

yi(x) = A cos([co> + 2nix]x) + B sin([c<j + 2nix]x) 
then satisfies the equation, 

-y'l = [u + 2mr] 2 y 1 , 
and at the endpoints of the edge we have 

y(0) = yi (0), y(l) = yi (l). 
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The endpoint derivatives are 

y'(0) = coB, y'(l) = uj\B cos(cj) - A sin(w)] 

and 

y [(0) = [uj + 2nrc}B, y[(l) = [u + 2nix] [B cos(u;) - A sin(cj)] , 

so 

, , , uj + 2nn ... ... lu + 2nn ... 

2/ = y'(0), y l = 2/(1). 

Thus yi satisfies the same interior vertex conditions that y does. More- 
over, the linear map taking y — > yi is one-to-one, as is the analogous map 
from E([u + 2nir} 2 ) to E(oo 2 ). 

□ 



2.2 Role of the discrete Laplacian 

Given a vertex v G Q, let n 1; . . . ,Ud eg ( v ) be the vertices adjacent to v. A 
discrete graph carries a number of linear operators acting on the vertex space 
H of functions / : V — ► C, including the adjacency operator 

deg(v) 

Af(v) = ]T /(«0. 
i=i 

and the degree operator 

Tf(v)=deg(v)f(v). 

Define the operator Ai by 

A 1 f(v) = f(v)-T- 1 Af(v). 
Ai is similar to the much studied Laplacian [H p. 3], 

The distinction of the cases A G {n 2 7r 2 } is related to the following fact. 

Lemma 2.2. Fix A G C, and consider the vector space of solutions of —y" = 
Xy on the interval [0, 1]. The linear function taking y(x) to (y(0),y(l)) is an 
isomorphism if and only if A ^ {n 2 7r 2 | n — 1, 2, 3, ... }. 
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Proof. For A ^ {n 2 7r 2 }, the formula 

y(x, A) = y(0) cos(a;x) + [y(l) - y(0) cos(a,)]^^ (2.4) 

shows that the map is surjective. On the other hand, if A G {n 2 7r 2 }, then 
y(0) = implies y(l) = 0, since y(x) = Bsin(nTrx). 

□ 

As an immediate consequence we have the following result for graphs. 

Lemma 2.3. Suppose the edges of Q have length 1, and A {n 2 ir 2 | n = 
1, 2, 3, . . . }. Let y : Q — > C be continuous, and satisfy —y" = \y on the edges. 
If y( v ) — at all vertices of Q , then y(x) = for all x G Q . 

Theorem 2.4. Suppose A ^ {n 2 ir 2 | n = 1, 2, 3, . . . } and y is an eigenfunc- 
tion for A^. If v has adjacent vertices u±, . . . , Ud eg (v), then 

deg(v) 

cos(u)y(v) = V 2/(ui). (2.5) 

deg(f) 



i=i 



Proof. In local coordinates identifying each ui with 0, (12.41) for yi on the 

edges Ci = (ui,v) gives 

• / \ / \ r / \ / \ / \ t ^ COS( Cl> J , , 

y {v) = -wsm w i/i iii + Wi(f) - yj(iii) cos(a;) , . (2.6) 

sin(w) 

Summing over i, the derivative condition at v then gives 

En \ i \ \ / \ Co> COs(c<j) v ■> 

y.(u) = -usin(w) > ^(ui) + . 7 r IjyAv) - yiiui) cos(o;)J. 

J J 2 

Using the continuity of ?/ at v and elementary manipulations gives (I2.5p . □ 

(12.51) is clearly an eigenvalue equation, with eigenvalue cos(a;), for the linear 
operator T _1 A acting on the space of (real or complex valued) functions on 
the vertex set. With this background established, we are ready to relate the 
spectra of A x and A^. 

Theorem 2.5. If X ^ {n 2 ir 2 | n = 0,1,2,...}, then A is an eigenvalue of 
Aqo if and only if 1 — cos(u;) = 1 — cos(vT) is an eigenvalue of Ai, with the 
same geometric multiplicity. 
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Proof. Since Ai = I — T~ 1 A, we may work with T~ 1 A. Suppose first that 
y(x, A) is an eigenfunction of Aoo satisfying the given vertex conditions. The- 
orem 12.41 shows that the (linear) evaluation map taking y : Q — > C to 
y : V — > C takes eigenfunctions to solutions of (12.51) . By Lemma [231 the 
kernel of this map is the zero function, so the map is injective. 
Suppose conversely that y : V — » C satisfies 

T _1 Ay(v) = \fi\ < 1. 

Pick A G cos _1 (/i). By Lemma 12.21 the function y : V — > C extends to a 
unique continuous function y(x, A) : £ — ► C satisfying — y" = Ay on each 
edge. 

In local coordinates identifying v with for each edge e; = (v, Ui) incident 
on v, this extended function satisfies (12.61) . Summing gives 

E„ s W r . 0/ \ 2/ m i \ ST~^ i \ W COs(co>) 

Vi{v) = ihiM [_sm {uj) ~ cos {uj)] 2^ Vi ^ Ui) + 2^) sinM 

i ^ ' i i ^ ' 

u; , s . , , , , cucosfu;) 



J^yi(wi) + deg(w)y(f 



sm(o;J sin(w) 
The vertex values satisfy (I2.5p . so 

En \ u ,uj cos(uj) 

Vi\ v ) = — de S( v ) cos{u)y{v) + deg{v)y{v) . = 0. 

sin(w) sin(w) 

Thus the extended functions are eigenfunctions of A ra satisfying the standard 
vertex conditions. Since the extension map is linear, and the kernel is the 
zero function, this map is also injective. 

□ 



2.3 Eigenfunctions at n 2 7r 2 

Now we turn to eigenvalues A G {n 2 ir 2 }. First recall [7J p. 7] that for both Ax 
and Aqo is an eigenvalue whose eigenspace is spanned by functions which 
are constant on connected components of Q. 

For n > 1 these eigenspaces for A^ also have a combinatorial interpreta- 
tion, closely related to the cycles in Q. If C is a cycle, and x is distance along 
the cycle starting at some selected vertex, then the function sin(2n7ra;) is an 
eigenfunction of Aqo. Similarly, if C is an even cycle, then sin(n7nr) gives a 
similar eigenfunction. We can make the following observation. 
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Lemma 2.6. Suppose Q is connected, andifj is an eigenfunction of A,^ with 
eigenvalue A = n 2 7r 2 for n > 1 . If ip vanishes at any vertex, then ip vanishes 
at all vertices. 

Proof. Suppose ip(v) = for some vertex v. On any edge incident on v, the 
eigenfunction is a linear combination 

ip(x) = Acos(mrx) + B sm(nnx), v ~ 0, 

and clearly A = 0. Since all edge lengths are 1, at all adjacent vertices w, 
we then have ip(w) = Bsm(mr) = 0. By continuity of ip and connectivity of 
the graph, ip vanishes at all vertices. □ 

The next result explores the existence of eigenfunctions vanishing at no 
vertices. 

Lemma 2.7. If X = (2nn) 2 , for n = 1,2,3,..., then Aoo has an eigenfunc- 
tion vanishing at no vertices. 

If A = (2n — l) 2 7r 2 , n = 1, 2, 3, . . . , then has an eigenfunction van- 
ishing at no vertices if and only if Q is bipartite. 

Proof. If A = (2n7r) 2 , then the desired eigenfunction is simply cos(2n7rx) in 
local coordinates on each edge. 

Suppose Q is bipartite, with the two classes of vertices labelled and 1. 
Pick local coordinates on each edge consistent with the vertex class labels, 
and define the eigenfunction to be cos([2n — l]7nr). 

Suppose conversely that for some A = (2n — l) 2 7r 2 , there is an eigenfunc- 
tion ip vanishing at no vertex. Label the vertices v according to the sign of 
ij)(v). In local coordinates for an edge, 

i/j(x) = acos([2n — 1]tcx) + bsin([2n — 1]ttx), a ^ 0. 

Then if w is a vertex adjacent to v we see that ip(w) = —ip(v), showing that 
vertices are only adjacent to vertices of opposite sign, and Q is bipartite. 

□ 

In addition to the vertex space mentioned above, an edge space may 
be constructed using the edges of a graph as a basis (we assume the field 
is K or C). The edge space has the cycle subspace Z (Q) generated by 
cycles, with dimension N E — N v + 1 [5, pp. 51-58] or [101 pp. 23-28]. Let 
E (n 2 7r 2 ) C E(n 2 7r 2 ) be those eigenfunctions of vanishing at the vertices. 
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Theorem 2.8. 

dimZ (£) = dim£ (4nV). 

Proof. It suffices to prove the result for a connected graph. Suppose Z (Q) 
has dimension M. Picking a spanning tree T for there is [5, p. 53] 
a basis of cycles Ci, . . . ,Cm such that each Cj contains an edge ej £ T, 
with ej not contained in any other Ci. Fix n G {1, 2, 3, ... } and construct 
eigenfunctions fj = sin(27rnx) on the edges of Cj, and on all other edges. 
Here x denotes distance along the cycle starting at some selected vertex. If 
a linear combination ^2 a ifi i s 0> then for x G ej 

M 

= ^2 a ifi( X ) = a jfj( X )' j = 1, • • • , M. 
t=l 

There are x G ej where fj(x) ^ 0, so a 3 - = and the functions fa are 
independent. This shows dimi? (4n 2 7r 2 ) > dimZ (^). 

Now suppose that tp G E (An 2 ir 2 ) . After subtracting a linear combination 
^2 aifi we may assume that if) vanishes on all edges not in the spanning tree 
T. For every boundary vertex v of T, if) vanishes identically on all but one 
edge of Q incident on v, and by the vertex conditions it then vanishes on all 
edges incident on v. Continuing away from the boundary of the spanning 
tree, we see that if) is the function. 

□ 




Figure 2.1: Bowtie graph 



The proof of Theorem [231 provides a combinatorial basis construction for 
A similar construction using even cycles will provide independent 
elements of E ((2n — l) 2 7r 2 ), but these may not form a complete set. Consider 
the bowtie graph in Figure 2.1, whose cycles have length 3. An eigenfunction 
if) with eigenvalue 7r 2 can be constructed by letting if>(x) = 2sin(7ra;) on the 
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middle edge. The function ip then continues as — sin(7r:r) on the adjacent 
four edges, and sin(7rx) on the remaining two edges. Notice that on each 
edge the function ip is an integer multiple of sin(7r:r). 

A generalization of the even cycles will be used for a combinatorial con- 
struction of E ((2n — l) 2 7r 2 ). Let Z\ denote the set of functions / : S — > C 
such that _ 

£/(e)=0, veV. 

Theorem 2.9. The linear map taking f G Z\ to g(x) defined by 

9e(x) = /(e) sin((2n - l)nx), 

is an isomorphism from Z\ onto E ((2n — l) 2 7r 2 ). The subspace Z\ has an 
integral basis. 

Proof. Since sin((2n— 1)ttx) = sin((2n— l)7r(l— x)), the edge orientation does 
not affect the definition, so g is well defined. Clearly g(x) satisfies the eigen- 
value equation and vanishes at each vertex. The condition Yle~v f( e ) = ® f° r 
all v G V gives the derivative condition, so g e E ((2n — l) 2 7r 2 ) . Moreover 
the map is one to one. 

Suppose ip(x) G E ((2n — 1) 2 tt 2 ). Then if) e (x) = a e sin((2n — l)irx) on 
each edge, and because ip satisfies the derivative conditions we have 

^a e = 0, v G V. 

Define /(e) = a e to get a linear map from E ((2n — l) 2 7r 2 ) to Z ± , which is 
also one to one. This establishes the isomorphism. 

To see that Z\ has an integral basis, let e n be a numbering of the edges 
of Q, and let x n = f(e n ). The set of functions Z\ is then given by the set of 
x±, . . . , xn £ satisfying the iVy equations 

^2 Xn = °- 

This is a system of linear homogeneous equations whose coefficient matrix 
consists of ones and zeros. Reduction by Gaussian elimination shows that 
the set of solutions has a rational basis, and so an integral basis. 

□ 
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3 Graph refinements 



Now we introduce the notion of graph refinement for graphs whose edge 
lengths are 1. Let the original combinatorial graph be denoted Qi, with the 
operator Ai acting on the vertex space Hi. For each integer N > 1 we will 
define a graph Qn with vertex space H^r and operator A^ : H^ ^» Hjv by 
subdividing each edge e£ft into N edges. Pick local coordinates identifying 
e with [0,1], and labelling = xq and 1 = xn- For n = 1,...,N — 1, 
introduce points x n = n/N in the local coordinates. The new graph Qn with 
/ = iVy + (N — 1) N £ vertices will have the vertex set consisting of the vertices 
of Qi, together with the new vertices x n , for each edge e of Q 1 . The vertex 
x n is adjacent to £ n _i and x n+ i for n = 1, . . . , N — 1. If v is a vertex in Qi, 
and the local coordinates for edges incident on v are chosen so v is identified 
with on each edge, then v is adjacent in Qn with the vertices x\ for each 
of the incident edges. An example is illustrated in Figure 3.1. 



Qi Q2 Qi 

Figure 3.1: Refinement of a graph 

An inner product on the vertex space H^v is defined by 

</' 9)n = ^ E deg(v)f(v)^), W = ^ deg(v) = 2NN £ . 

Using the above identifications of edges of Q N with subintervals of [0,1], 
we may also construct the corresponding continuous graph having edges of 
length 1/N. These continuous graphs may all be identified with the contin- 
uous graph Goo corresponding to Q\. In particular they will share the inner 
product 

(f,g)oo = tt / fa, 
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where N £ is the number of edges in Q\. As a continuous graph, the standard 
vertex conditions hold at all vertices of Qn- Although appearing different in 
definition, the continuous graph Laplacian A^ has not changed as we pass 
from Qi to Q N . 

Let Atv denote the following operator on the vertex space of Qn, 

A N = N 2 (I -T^A). 

The eigenspaces of these operators with eigenvalue A will be denoted E N (\), 
while £"00 (A) will denote an eigenspace for A^. By ( I2.3P and Theorem 12. 5\ 
the mapping iV 2 (l — cos(V\/N)) carries eigenvalues of Aoo to eigenvalues 
of Atv if X/N 2 {n 2 rr 2 }. That is, the eigenvalues of the normalized adja- 
cency operator are cos(\/A/iV) whenever A ^ {(Nun) 2 } is an eigenvalue for 
Aoo. Notice the eigenvalues of An are approximately A/2 when A is small 
compared to N. 

Lemma 3.1. For 1 < N < 00 the operator An is self adjoint on the vertex 
space with the inner product (f,g)N- 

Proof. It suffices to check the normalized adjacency operator T~ X A. □ 

Let E p (n 2 ir 2 ) denote the subspace spanned by eigenf unctions of hav- 
ing the form cos(nirx) on each edge (so not vanishing at the vertices). Let 
Sat C L 2 (^oo) denote the subspace 

§jv = span{E p (A^ 2 7r 2 ), J E oo (A),0 < A < iVV}. 

Here is preliminary result describing the restriction of eigenf unctions of A^ 
to the vertices of Qn- 

Proposition 3.2. The restriction map Rn : §at — > Mn is an bisection. For 
< A < N 2 it 2 this map takes distinct orthogonal eigenspaces £ , o(A) of A^ 
onto distinct orthogonal eigenspaces E N (N 2 (l-cos(V\/N))) of A N , and R n 
takes E p (N 2 n 2 ) onto E N (2N 2 ). 

Proof. By Theorem I2.5I and (12. 3ft . Rn is a bijection from the eigenspace 
-^00 (A) of Aqo to the eigenspace E N (N 2 (1 - cos( v / A/iV))) of A N . The 
eigenspaces for A^ and An are just the functions which are constant on the 
connected components of the respective graphs, so Rn is a bijection from 
£"00(0) to £jv(0). From the proof of Lemma [27fl we also see that Rn is a 
bijection from E p (N 2 tc 2 ) to E N (2N 2 ). 
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Since and A N are self adjoint on their respective function spaces, 
and since cos(t) is strictly decreasing on (0,7r), distinct eigenspaces of A^, 
which are orthogonal, are mapped to distinct eigenspaces of A at, which are 
orthogonal and span H^. □ 

The restriction map Rn '■ §jv - ► Hjv also has noteworthy features on the 
individual eigenspaces E^X) of A^. Before stating the results, we start with 
a simple observation relating the EItv inner product and the sums appearing in 
the trapezoidal rule for integration. Recalling that W = J2 V deg(v) = 2NN £ , 
the inner product for Qn satisfies the identity 

W(f,g) N = J2 de 9(v)f(v)W) 

V 

N-l 

= ^2[fe(x )g e (x ) + fe(x N )g e (x N ) + 2 ^ f e (x n )g e (x n )}. 

eeGi n=l 

These last sums are just the trapezoidal rule sums used for integrals over 
the edges of Goo. With this motivation, if the continuous linear functional 
T N : C[0, 1] -> C is defined by 

1 

T N (f) = 2]y[/(so) + f(x N ) + 2 /Ml, 

n=l 

then 

(R N f,R N g) N = j^J2 T *(f^)- (3- 1 ) 
The following identities are also useful. First 



L 



1 e 2iu,x = exp(2^) - 1 = eluJ sin{u) 



Then, for < u < N, a geometric series computation gives 

N-l 

2NT N (exp(icux)) = exp(icu) — 1 + 2 ^ exp(iun/N) 

n=0 

= exp (iu) - 1 + 2- ;. v .; n = 1 - expM)- P . 'J. , 
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or 



e™, Mo{z) = Mo(-z) = acot(z). (3.2) 

These identities will be useful for comparing the inner products on L 2 ({?oo) 
and H.N- The cases A = k 2 7r 2 with < k < N are considered first. 

Theorem 3.3. Suppose f,g G E^X), with < A < N 2 n 2 . If A = k 2 n 2 for 
an integer k with < k < N, then 

(f^g)™ = (Rn/i Rn9)n- (3.3) 

Iff,ge E p (N 2 tt 2 ), then 

(f,g)oo = \(RNf,R N g)N- (3.4) 

Proof. Suppose /, g G E^X), with < A < N 2 ti 2 . On each edge e we have 
f e = a e e^ x + f3 e e- l ^ x , g e = le e^ x + ^e"^, cu 2 = A, (3.5) 

/e£ = «e7l + + ^Ae 2 ^ + PeTeP'^ ■ 

If A = A; 2 7r 2 for an integer k with < k < N, then 



T N {e 2 ™ kx ) = [ e 2nikx = 0. 
Jo 



Using (13.11) and T N (1) = 1, it follows that 

If i 

(f,g)oc = ^ fg = ^y>ei£ + /3A] 

= --7- ^ T N (a e % + /3 e t? e ) = (R N f,R N g) N , 

c e 

which is (13.31) . Similar calculations handle the case A = 0. 

Suppose f,g G E p (N 2 tt 2 ), a one dimensional space. For / e = g e 
cos(iV7rx), we have 



1 



cos 2 (A%x) = 1/2, 
15 



but 



giving (IcTIj) . 



□ 



For general co> a variation of a classical trapezoidal rule estimate appears. 

Theorem 3.4. Suppose f,g e E^X), with < A < iV 2 7r 2 . TTien /orw > 
u; > 0, 

(i2jv/,^>jv -</,«/>! =0(|u;| 4 /iV 4 )||/|| \\g\\. (3.6) 

Proof. The argument starts with a variation of standard error estimates [U p. 
285], [p. 358-369] [6] for the trapezoidal rule for integration. If = exp(iujx) 
then 



0'(1) - 0'(O) 



12iV 2 12iV 2 
Putting this together with (I3.2p yields 



-a; 



12iV 2 



e —l N {e ) 



12iV 2 U 2iV' 



with 



Mx(2) = 1 - scotte) - 



A Taylor expansion gives 

M 1 (z)=0(z A ), 

and the function M\{z) = M\{—z) is analytic for \z\ < tt. 

On each edge we have the representation (13.51) . Since / and g are in the 
domain of Aoo, the function ip = f~g satisfies the vertex conditions (12.11) and 
( 12T2D . Thus 

Et^(i) - ^(o)] = -EE a^oo = °> 

e 1>GV e~i; 

and 



E 



N £ [(f,g)-(R N f,R N g) N ]=J2\ [ ^-T N {^ e 

^(i)-^(o)- 



12iV 2 



E + 



7e 



16 



Mil 



N J 



1 

o ~~ •/() 

2 i 1/3 12 



Using 2|a e ||/? e | < \a e \ 2 + |/3 e | , we find 



\f e \ 2 = \a e \ 2 + \(3 e \ 2 + af3e 



a; 
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>(1- 

This gives the desired result 



sm w 



a; 



K| 2 + |& 



< 



(R N f, R N 9)n - (f,g) 

f 1 \fe\ 2 ) l/2 I \9ef) 1/2 



< 



-'2(/ l/D 1/2 (/ IjI 



2U/2 



□ 



4 The FFT 

In this section we turn to an algorithm for the efficient computation of the 
Fourier transform and its inverse transform on the spaces M.^. 







7T 



2tt 



3tt 



47T 



57T 



Figure 4.1: Square root of the spectrum of Aoo. 

Recall some of the structures associated with the spectrum of A^. The 
square roots V%- of the positive eigenvalues of and the eigenspaces E(Xk) 
exhibit the '27r -periodicity' discussed in Proposition 12. 11 and illustrated in 
Figure 4.1. Let < co> ,i < u; ,2 < ■ ■ ■ < 2n denote the square roots of 
the K distinct eigenvalues of A^ with < < (2n) 2 . For each ojq^ 
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pick an orthonormal basis ^(uj ^) for E(ujQ k ), with basis functions ^(ufyfc), 
j = l,...,dim£?(a;g fc ). 
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Figure 4.2: Eigenfunction maps for Aoo. 

For m — 1, 2, 3, . . . , let uj m ^ = ^o,fc + 2mTr. From each basis ^(uo k) we 
may produce bases *f?(u) m ,k) f° r by the method of Proposition 12.11 

Recall that on each edge e a function tyj e ty(wo t k) has the form 

tyj = A e exp(iu> Qt kx) + B e exp(—iuio t kx), < x < 1. 

The set of functions *f?j(u mt k) whose values on e are 

= A e exp(i[u 0yk + 2m7[}x) + B e exp(-i[(j ,k + 2m-n]x), (4.1) 

then give a basis ^{oJ m ,k) f° r ^([^o,* + 2m7r] 2 ). This basis need not be 
orthonormal. Figure 4.2 indicates the relations among the bases ^>(uj m ^)- 

Proposition 4.1. If £ { 7r ?27r}, then the functions ^ j{uj m ^) are or- 
thonormal, as are the functions R^^> j(oJ m ,k) £ ^-n f or < uj m ^ < Ntt. 

Proof. Suppose ^/i(u mk ) e f(w m]k ) is represented on the edge e by 
^i(^m,k) — C e exp(i[ui 0t k + 2mii}x) + D e exp(—i[u Qk + 2rrni}x), < x < 1, 
and (14. ID describes tyj. With uj ^ = vr,27r, 



/ exp(±2i[u; 0]fc + 2rriTc]x) = 0, 



so 



; f Vjiu^MVm*) = Y\[^P e + B e D e ] = V / ^(c^A^K 

; : Jo 



e " " e e 



The extension to Hw follows from Theorem I3J31 □ 
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Let I = dim_E(0), Ik = dim.E(ujQ k ), and let 

be the restriction of ^j{uJ m ,k) to the vertices of Qn- The bases <&{uJk,m) will 
be used for efficient Fourier transform algorithms [U p. 182] [HI p. 383]. 
Define the discrete Fourier transforms (DFT) Tn : H^v — > C /() @ m ,k < & Ik f° r 
the continuous graph by taking the inner product of a function / e M. N 
with the bases ${oj m ^)-, 

FnU) = {(f, ^m,k)) N }, u m , k < Nn. (4.2) 

The condition uj m ^ < iW amounts to m = 0, . . . , iV/2 — 1 when N is even, 
which will hold in the cases of interest below. Except for the cases noted 
in Proposition 14.11 the functions of $(u; mj fc) may not be orthonormal, so the 
Fourier transform will not be an isometry from $(o; TO> fe) to C /fc if the usual 
inner product is used on the range. We consider a modified inner product. 

Suppose B{ijj m k ) = {pij) is a matrix taking the basis $(u; m) &) to an or- 
thonormal basis 

Vi = /Mj®A u m,k) (4.3) 

3 

for the same eigenspace in Ef/v. (Such a matrix may be obtained by the Gram- 
Schmidt process.) For / in the span of $(u; TO> fc), the map / — > {(/, tj^n} is 
an isometry to C /fc with the usual inner product X • Y = ^ x jVj- Expressing 
this in the original basis, 

imr = E \(f>viU 2 = e (/• >>/!•,>. v 2 = E i EW. ^)m 2 

i i j i j 

= BX.BX, X = ({f, ^) N , ...,(/, $ Ik ) N ) T . 

For / = J^. Cj$j(o; mj fc) in the span of $(u; mj fc), the coefficients Cj may be 
recovered from the DFT values (/, $j)tv. Starting from 

/ = E ^' fli = (f'^N, 
i 

and using (14. 3[) we find 

i j 2 I j i 
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l 3 

or 

C = X)(S'B) iI (/ J $ i V (4.4) 

3 

Calling BX • BX the B(u mt k) inner product on C /ft , and noting 3>(o> m) fc) 
is a basis for the \i m ^ eigenspace of A N , we obtain the first part of the next 
result. 

Theorem 4.2. The Fourier transform T N : M N — > C /o © mi fc C /fc satisfies 

•7"iv(Ajv/) — {^m^ifjm^fi 

and is an isometry if /ias t/te B(u m k) inner product. 

If N is a power of 2, then ^tv(/) can be computed in time 0(Nlog 2 (N)). 

Proof. The inner products of f!4.2[) may be split into trapezoidal sums for the 
edges as in (13. ip . 

£ e 

Using the representation (14.11) we have 

T N (f e [A e exp(i[uj 0tk + 2rmr]x) + B e exp(-i[u 0)k + 2ra7r]a;)]) 

= A e T N (exp(—iuj ^x) f e exp(—2mnix) + B e T N (exj>(iuJo j kX)f e exp(2m7Tix)). 

If N is a power of 2, e and k are fixed, and m = 0, . . . , y — 1, then 
the sequence T^(exp(±iujQ^x)f e (x) exp(±2rrnrix)) differs trivially from the 
conventional discrete Fourier transform of the sampled data f e (x n ), which 
may be computed pQ p. 182] [6, p. 383] in time 0(N\og 2 (N)). Since the 
number of edges e and indices k is fixed, we obtain the result. 

□ 

It remains to consider the inverse transform. Since the matrix B*B of 
(14.41) depends only on the graph and functions <&j (a> m ,&) , these matrices are 
assumed to be precomputed. Suppose the DFT sequence f!4.2[) of a function 
/ is given. Since the computations in (14.41) are limited to sequence blocks 
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C /fc , the coefficients Cj(uj m ^) for all may be computed in time O(N). It 
remains to compute 

m,k 

Fixing the edge e, the primitive frequency index k, and the basis index j, 
and taking advantage of (14.1 1) , the remaining sum has the form 

f e ,k,j( x n) = A e exp(iu 0tk n/N) Cj(u; mjfe ) exp(i2rmm/N) 

m 

+B e exp(— iujQ^n/N) Cj{uj m ^) exp(— i2m7rn/N). 

m 

These sums may be computed with a conventional FFT. Since the number 
of edges e and primitive frequency indices k are constant, and the range of 
basis indices j is bounded independent of m, the next result is obtained. 

Theorem 4.3. If N is a power of 2, then the inverse DFT can be computed 
in time 0(Nlog 2 (N)). 



5 A family of examples 

A relatively simple family of examples is provided by the complete bipartite 
graphs K(m, 2) on m, 2 vertices. The m = 4 case is portrayed in Figure 
5.1. Discussion of these graphs is facilitated if the vertices are colored, the 2 
vertices of degree m being red and the remaining m vertices blue. 

Suppose v is a vertex, and {ui} are the vertices adjacent to v. Then the 
operator A x on the m + 2 dimensional vertex space is given by 

deg(v) 

Eigenvalues \l for Ai are 0,1,2, with 1 having multiplicity m. The eigen- 
functions ip for /x = 0, 2 are 

yUo = 0, lp (v) = 1, 
/Wi = 2, = {_\ ^bte} 
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Figure 5.1: The graph K4 2 
Next, consider the eigenf unctions for // = 1. A function f(v) will satisfy 

) = f(v) if and only if Ylt^i"^ f( u i) — f° r every vertex v. Number 
the red vertices r\ and r 2 , and the two blue vertices bi, . . . ,b m . For j = 
1, . . . , m — 1, independent eigenfunctions are 

1, f = 6j 

^» = {-1, u = 6j+i,}. 
0, otherwise 

There is one additional independent eigenfunction 

1, v — Ti, 
ip m (v) = u = r 2 , }. 

0, otherwise 

To establish the independence of ^1, ... , note that if 

m 

then c m = by evaluation at r l5 and evaluation at the points bj leads to 
equations 

Ci = 0, Ci = c 2 , . . . , c m _i = c m . 

Moving to the continuous graph, identify each edge with [0, 1] so the blue 
vertex is identified with 0. Consider the eigenvalues A and corresponding 
eigenfunctions ^ for Aoo, focusing on y/\ e [0, 2ir\. Of course we have 
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Ao = with the constant eigenfunction \l/(0) = 1. From fij = 1 we obtain 
eigenvalues 

= cos -1 (l - fy) = cos _1 (0) = tt/2, 3tt/2. 

That is, A = (tt/2) 2 and A = (37r/2) 2 are eigenvalues, each with multiplicity 
m. The corresponding eigenfunctions tyj(jr/2) and ^/j(3ir/2) are obtained by 
extrapolation as in Lemma [2.21 from the vertex values of ipj. 

This graph [5, p. 53] has m — 1 independent cycles, all of which are 
even. One basis of cycles Cj is given by the vertex sequence n, r%, foj+i, r\. 
Each Ci supports eigenfunctions ^fi(ir) and \E^(27r) which are respectively 
± sin(7rx) and ± sin(27rx) on the edges of the cycle. 

Finally, there are two additional independent eigenfunctions, \& m (7r) hav- 
ing values cos(7ra) and ^ m {2n) with values cos(27rx) on the edges of Q. This 
gives a total of 

2N £ = 4m = 2m + 2{m - 1) + 2 

independent eigenfunctions for Aoo with eigenvalues < A < 2n. 

This graph has an obvious involution in which r± and r 2 are interchanged, 
and we may split the eigenspaces into even and odd subspaces with respect to 
the involution. Notice that the eigenfunctions \& m (7r), \& OT (27r), \fri(ir/2) and 
*i(37r/2) for i = 1, . . . , m - 1 are even while ty m (ir/2), ^ m (37r/2), and 
^j(27r) for i = 1, . . . , m — 1 are odd. By restricting to the even, respectively 
odd, subspaces we may treat the case of a star graph with m edges and 
boundary vertices with the conditions f' e (0) = 0, respectively / e (0) = 0. 
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